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Large-scale quantum computers must be built upon quantum bits that are 
both highly coherent and locally controllable. We demonstrate the quantum 
control of the electron and the nuclear spin of a single atom in silicon, using 
a continuous microwave magnetic field together with nanoscale electrostatic 
gates. The qubits are tuned into resonance with the microwave field by a local 
change in electric field, which induces a Stark shift of the qubit energies. This 
method, known as A-gate control, preserves the excellent coherence times and 
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gate fidelities of isolated spins, and can be extended to arbitrarily many qubits 
without requiring multiple microwave sources. 

The construction of a large-scale quantum computer is among the most exciting scientific 
challenges of our time. For this purpose, the ability to exploit the advanced fabrication meth¬ 
ods developed in semiconductor nanoelectronics {1,2) would be highly beneficial. Succeeding 
at this project will depend upon the ability to create quantum bits (qubits) which are, at the 
same time, highly coherent and easy to control and couple to each other. However, quantum 
coherence and ease of control are often antithetic requirements. Several types of semiconductor 
spin qubits have been designed to be operated using only electric fields (3, 4,5, 6). This allows 
control at the nanometer scale via small electrodes, but can lead to unwanted decoherence from 
charge and gate noise (7, 8). Other spin qubits, such as the electrons and nuclei of donors 
in silicon (9, 10) or nitrogen-vacancy (NV) centers in diamond (11,12), exhibit outstanding co¬ 
herence times thanks to their lack of sensitivity to electrical noise and the reduced nuclear spin 
fluctuations in their vicinity. However their coherent operation requires high-frequency oscillat¬ 
ing magnetic fields, which are difficult to confine to nanometer scales. A solution to reconcile 
long coherence with local electrical control was already described in the visionary proposal 
of Kane (13) for a donor-based quantum computer in silicon. The resonance frequency of a 

donor spin depends on the applied magnetic field Bq, as well as the electron-nuclear hy- 
perfine coupling A. The latter can be locally tuned with an electrostatic gate, known as the 
“A-gate”, which deforms the wave function of the donor-bound electron and modifies the prob¬ 
ability density at the nucleus (14). Kane envisaged a multi-qubit quantum computer where a 
global, always-on microwave magnetic field is by default off-resonance with the qubits. Quan¬ 
tum operations are controlled by locally modifying A in order to bring the desired qubits into 
resonance with the global microwave field (15). This eminently scalable proposal has been 
further developed to include quantum error correction (16). Here, we present the experimental 
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demonstration of local and coherent electrical control of both the electron and the nuclear spin 
of a single donor in silicon, representing the first realization of Kane’s A-gate on a single 
qubit. 

The donor in silicon constitutes a two-qubit system, where both the electron (indicated 
with 14,) or It)) and the nuclear (|JJ,) or jfl')) spin states can be coherently controlled by a mag¬ 
netic field Bi oscillating at specific electron spin resonance (ESR) and nuclear magnetic res¬ 
onance (NMR) frequencies. We fabricated a device that comprises a single ^^P donor in an 
isotopically purified ^®Si epilayer (17), implanted (18) next to the island of a single-electron- 
transistor (SET). The SET is formed under an 8 nm thick Si02 layer by biasing a set of elec¬ 
trostatic gates (yellow in Fig. lA). The distance between donor and SET island is ~ 20(5) nm, 
resulting in a tunnel coupling of order 10 kHz. The device is cooled by a dilution refrigerator 
(electron temperature Tei 100 mK), and subject to a static magnetic field Bq — 1.55 T applied 
along the [110] Si crystal axis. Due to the Zeeman effect the electrochemical potential /r of the 
donor electron depends on its spin state, with > /ij,. Another set of gates (pink in Fig. lA) 
is used to tune the electrochemical potentials of donor and SET island (//set) to the readout 
position (“Read/Init.” in Fig. 2A), where > //set > l^i and only the |t) state can tunnel 
out of the donor. The positive donor charge left behind shifts the SET bias point and causes a 
current to flow, until a jj,) electron tunnels back onto the donor. This spin-dependent tunneling 
mechanism is, therefore, used to achieve single-shot electron spin readout (19,20), as well as 
14,) initialization. For coherent spin control, an oscillating magnetic field Bi is delivered to the 
donor by an on-chip, broadband transmission line terminating in a short-circuited nanoscale 
antenna (27) (blue in Fig. lA). While manipulating the electron spin state, the gates are tuned 
such that < uset, to ensure that the electron cannot escape the donor (“Pulse ESR/NMR” 
and “A-Gate Control’’ positions in Fig. 2A). 

In the present experiment, the same gates used to tune the donor to the readout position are 
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Figure 1: Electric field dependence of electron and nuclear energy states. A, False-colored 
scanning electron microscope image of a device similar to the one used in the experiment. Blue: 
microwave (MW) antenna; Yellow: gates used to induce the single-electron transistor (SET) 
charge sensor under the Si02 insulator; Pink: “A-gates”, comprising gates labeled Donor Fast 
(DF), Donor Slow (DS) and Top Gate AC (TGAC). These gates are used to tune the potential and 
electric field at the donor location. B, Electron wavefunction of a donor under an electrostatic 
gate. A positive voltage applied to the gate attracts the electron towards the Si-Si02 interface. 
Eor illustration purposes, the wavefunction distortion is largely exaggerated as compared to the 
actual effect taking place in the experiment. C, Energy level diagram of the neutral e“-^^P 
system. Gate-controlled distortion of the electron wavefunction modifies A and je, shifting the 
electron spin resonance (ESR) i/gi and r'e 2 , and the nuclear magnetic resonance (NMR) i^ni and 
r '„2 transition frequencies. 


also used to apply an electric field shift, which in turn modifies the spin transition frequencies 
via the Stark effect. Therefore, we label them “A-gates” henceforth. In Pig. IB we illustrate 
the effect of the electric field on the donor electron. A positive bias on a gate located above 
the donor pulls the electron wave function towards the Si/Si 02 interface and away from the 
nucleus (J4). This modifies both the hyperfine coupling A (~ 117.53 MHz in bulk, and 96.9 
MHz for this device) (14,22) and the electron gyromagnetic ratio 7 e (~ 27.97 GHz/T in bulk) 
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Figure 2: Local electrical control of the ESR transition frequencies. A, Schematic of the 
charge stability diagram for this device. The thick solid lines represent the Coulomb peaks of 
the SET, while the dashed line indicates the ionization/neutralization of the donor. B, Measured 
shift in ESR frequencies i^ei, 2 (VA) as a function of the A-gate voltage Va- Accurate values of 
i^ei, 2 (VA) are obtained by coherent Ramsey experiments (see text). The change in local electric 
field is obtained from finite-element electrostatic modeling for the specific device geometry and 
donor location (see supplementary section S6 for details). 


(23). The ESR frequencies 1 / 61,2 depend on both these parameters, and on the state of the 
nuclear spin. In the limit 7ei?o ^ we find z/gi == 7e^o — ^/2 (for nuclear spin |JJ-)) and 
i^e 2 = 7e-Bo + ^/2 (for 111)), as shown in Fig. 1C. Identifying whether the instantaneous ESR 
frequency is i/gi or i^e 2 constitutes a single-shot, quantum nondemolition nuclear spin readout 
(JO). The NMR transitions can be gate-tuned as well, albeit only via the Stark-shift of A, 
since i^ni = A/2 + 'ynBo (for electron spin j).)) and 12^2 = A/2 — 'jnBo (for |t)), where 
7 „ = 17.23 MHz/T is the nuclear gyromagnetic ratio. 
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During pulse-ESR/NMR experiments we normally operate the electrostatic gates in a com¬ 
pensated manner, to keep /Uset constant while shifting with respect to it (79) (“Pulse 
ESR/NMR” position in Fig. 2A). This, however, results in a limited variation of the electric 
field at the donor site. To induce a significant Stark shift we adopted an uncompensated gat¬ 
ing scheme, where both and //set are drastically lowered (“A-Gate Control” position in 
Fig. 2A, see supplementary section S4 for details on the exact gate configuration). The re¬ 
sulting change in electric field can be calculated with a finite-element Poisson equation solver 
(TCAD) for the specific device geometry and the triangulated donor location (see supplemen¬ 
tary sections S5 and S6 for details) (24). 

We demonstrate the gate-induced shift of the ESR frequencies by performing, with conven¬ 
tional pulse-ESR, a series of Ramsey experiments on the electron spin at different values of Va 
(Fig. 2B), where we define = 0 as the readout position. The spin is rotated from |).) to the 
xy-plane by a 7r/2 pulse at frequency i^mw> then left to freely precess for a time r, then rotated 
by 7r/2 again. The accumulation of a phase shift between the spin precession at Ue{VA) and the 
MW reference clock at z^mw gives rise to oscillations in the probability of finding the electron 
It) at the end of the sequence. The frequency of the Ramsey fringes gives a very accurate value 
for ^^{Va) — I'Mw- 

Both z/ei and z/e 2 shift to lower frequencies upon increasing Va (Fig. 2B). This indicates 
that a significant Stark shift of % (i.e. the electron p-factor) is taking place, in addition to the 
A-shift. Linear fits to z^ei, 2 (VA) yield slopes di^ei/dVA = —2.27(6) MHzW and dPe 2 ldVA = 
—1.36(3) MHz/V. Using the expressions given in Fig. 2B we extract the tuning parameters 
aA = dA/dVA = dUe 2 /dVA — dUei/dVA = 0.91(7) MHzW and = d'yeBo/dVA = 

{due 2 /dVA + duei/dVA)/2 = -1.81(5) MHz/V, with = d^iJdVA = -1.17(3) MHz/V/T 
at Bq = 1.55 T. The positive value of a a indicates that increasing Va leads to an increase in 
the electron probability density at the nucleus. This is because a strong electric field is already 
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present for the purpose of forming the SET, whereas increasing Va causes an additional electric 
field in the opposite direction, thus an overall reduction of the hyperfine Stark shift (see finite 
element simulations in supplementary section S6). The absolute value of the hyperfine coupling 
and its tunability is in agreement with atomistic tight-binding simulations with NEMO-3D (25, 
26) for the range of donor positions, electric fields, and strain expected for this device (see 
supplementary section S7). 

An A-gate voltage Va ~ 300 meV results in a frequency tuning range = 400 — 

700 kHz. Thanks to the use of an isotopically enriched ^®Si epilayer (77), this is a factor > 200 
larger than the intrinsic linewidth g ^^2 of the ESR transitions (22). Eor the same 

Va, the NMR frequencies can be shifted by = 125 — 150 kHz by the Stark shift of A 

alone, which is a factor r\j 250 higher than the intrinsic linewidth = 0.5 kHz (22). 

Having calibrated the voltage-controlled qubit frequency shifts, we demonstrate how to use 
A-gate pulses to perform coherent control of the qubit states around the Bloch sphere in the 
presence of a continuous-wave (CW) oscillating magnetic field Bi. We demonstrate this on both 
the electron and the nuclear spin qubits. Eig. 3A is a schematic of the electrically controlled 
Rabi sequence. We pulse lA to V7, the voltage needed to tune the spin transition in resonance 
with the MW or RE source, for a time tp, to coherently drive the spin around the A-axis of the 
Bloch sphere, defined in the reference frame rotating at i^mw (for ESR) or i/rf (for NMR). The 
coherent Rabi oscillations (Fig. 3B,C) are used to calibrate the duration of the control pulses for 
any desired rotation angle. 

An electrically-controlled Ramsey experiment (Fig. 3D) is obtained by tuning the spin tran¬ 
sition frequencies into resonance for the duration of a 7r/2-rotation, then moving them to a 
detuned value i^det for a wait time r to accumulate a phase shift 4*^ = t(z/(K-) — z/det) with 
respect to the ly(Vr) reference frame, and finally tuning them back into resonance for the second 
7r/2-rotation (Fig. 3E,F). Qubit rotations around the E-axis can be achieved by accumulating 
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Figure 3: Electrically-controlled qubit control and coherence measurements. A, Schematic 
of the sequence used to measure EC Rabi oscillations. B, C, EC Rabi oscillations measured on 
the electron and nucleus, respectively. D, Schematic of the sequence used to measure EC 
Ramsey oscillations. E, F, EC Ramsey oscillations measured on the electron and nucleus. 
G, Schematic of the sequence used to measure EC coherence times. H, I, EC Hahn echo decay 
for the ^^P electron and nucleus. J, K, Extended spin coherence times T 2 for CPMG dynamical 
decoupling sequences on the ^^P electron and nucleus. 


an additional = 7r/2 before bringing the spin into resonance. This allows for full two- 
axis (X and Y) control of the qubits’ states on the Bloch sphere. Alternatively, the phase 
accumulated while off-resonance can be used to produce a controllable Z-gate. 

We measure the qubits’ coherence times T 2 by performing Hahn echo and Carr-Purcell- 
Meiboom-Gill (CPMG) dynamical decoupling sequences (Pig. 3G-K). Por these measurements 
the wait times r are chosen such that they always result in a multiple of 27r phase shift. The 
Hahn echo yields = 0.97(6) ms and = 5.1(3) ms for the donor electron and nuclear 
spins, respectively (Pig. 3H,I). A CPMG sequence (see schematic in Pig. 3G for details) fur- 
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Figure 4: Electrically controlled gate fidelities. A, B EC randomized benchmarking per¬ 
formed on the electron and nucleus, respectively. Shaded circles are the results of individ¬ 
ual measurements (i.e. individual random sequences of gate operations), while the solid circles 
show the average state survival probability of all random sequences with the same number of 
gate operations. 


ther decouples the qubits from low-frequency noise, and extends the coherence times to up 
to = 10.0(8) ms and = 98(9) ms by applying 32 and 256 yr-pulses, respec¬ 

tively (Fig. 3J,K). The coherence times of the electrically-controlled electron qubit, measured 
in milliseconds, can be fitted by = 1.03A^° ®^. Conventional pulse-ESR operation 

(black diamonds in Fig. 3J) gives nearly identical values, = 0.93A^°'^° (see also 

supplementary section S10 for more data). The electrical control method preserves the excel¬ 
lent spin coherence of the qubits, because no additional coupling between the spin states and 
the environment has been introduced. 

We apply the two-axis control to construct the complete set of 1-qubit Clifford gates, and 
conduct randomized benchmarking (27, 28) experiments (Fig. 4A,B) to quantify the average 
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gate fidelities using electrically-controlled qubit manipulation (Fig. 4). The method relies upon 
the measurement of the probability Vi^N) of arriving at the correct final qubit state after per¬ 
forming a long sequence of randomly-chosen quantum gates. In the presence of gate errors, 
the probability decays according to the number of operations N. We then extract the average 
fidelity for a single Clifford gate by fitting the equation 

P(A^) = M(2jr^ - 1)^ + 0.5 (1) 

to the data {28). M is a free parameter, and depends on the initialization and readout fidelity. We 
obtain = 99.0(1) % and = 99.3(1) % for the electron and the nuclear spin, respectively. 
Since each Clifford gate is composed on average of 1.875 individual gate operations (29, 28), 
we can also quote average single-gate fidelities of = 99.4(1) % for the electron and 

jz^mgie _ 99 5 ( 1 ) % for the nucleus. These gate fidelities are comparable to those obtained by 
pulse-ESR/NMR randomized benchmarking for similar microwave powers {28), and are mostly 
limited by the ratio of gate time to coherence time. 

The speed of an electrically-controlled gate operation is inherently limited by the voltage 
tunability of the resonance frequencies. The excitation profile of the CW field (Fig. 3A) must 
be narrow enough to leave the qubits unperturbed while off-resonance, imposing the condition 
7 i?i <C Az/. For linearly oscillating Bi, {'^Bi)~^ is the duration of a 7r-rotation. Therefore, 
Ai/ < 1 MHz in this device requires gate times > 10 /iS. 

The qubit control method demonstrated here is applicable to any resonantly-driven qubit 
where the resonance frequency Vr can be quickly and locally controlled by an electric field, and 
shifted by much more than the resonance linewidth. A wide variety of qubits can potentially 
fulfil this condition. For spins in diamond {30) and silicon carbide {31) Vr can be tuned by 
modifying the crystal field parameters in the spin Hamiltonian. Magnetic molecules can have 
tunable Vr through a hyperfine Stark effect {32) similar to the one shown here. Several types 
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of semiconductor quantum dot qubits exhibit tunable electron spin ^(-factor (33, 34) through 
spin-orbit coupling effects, or tunable splitting through the interplay of valley-orbit and tunnel 
couplings (6). 

Due to the high cost of vector microwave signal generators, it seems implausible that future 
multi-qubit experiments will resort to a dedicated source for each qubit. Time- and frequency¬ 
multiplexing qubit control is possible in proof-of-principle experiments, but is impractical in 
large fault-tolerant quantum processors. In most error correction schemes, fault tolerance is 
only guaranteed if all qubits can be operated simultaneously at any time. The method of qubit 
control demonstrated here fulfills all the practical requirements for a large-scale quantum com¬ 
puter, since control gates can be applied simultaneously to arbitrarily many qubits, while requir¬ 
ing only one CW microwave source together with inexpensive multi-channel baseband pulse 
generators. 
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Supplementary Information 


51 Device Fabrication 

The device was fabricated on a 0.9 /im thick epilayer of isotopically purified ^®Si, grown on 
top of a 500 fim thick “^*Si wafer. The ^®Si has been depleted to 800 ppm in the enriched 
^®Si epilayer. Single-atom qubits were selected out of a small group of donors implanted in a 
region adjacent to the Single-Electron-Transistor (SET). In this device, molecular ions were 
implanted at 20 keV energy in a 100 x 100 nm^ window. All other nanofabrication processes 
were identical to those described in detail in Ref. (7), except for a slight modification in the gate 
layout to bring the qubits closer to the microwave antenna and provide an expected factor 3 x 
improvement in Bi (see Fig. S2 for schematic of the gate layout). 

52 Experimental Setup 

The sample was mounted on a high-frequency printed circuit board in a copper enclosure, ther¬ 
mally anchored to the cold finger of an Oxford Kelvinox 100 dilution refrigerator with a base 
temperature Tbath = 20 mK. The sample was placed in the center of a wide-bore superconduct¬ 
ing magnet, oriented so that the Bq field was applied along the [110] plane of the Si substrate, 
and perpendicular to the short-circuit termination of the MW antenna. The magnet was operated 
in persistent mode while also feeding the nominal current through the external leads. We found 
that removing the supply current while in persistent mode led to a very significant magnetic 
field and ESR frequency drift, unacceptable given the intrinsic sharpness of the resonance lines 
of our qubit. Conversely, opening the persistent mode switch led to noticeable deterioration of 
the spin coherence, most visible as a shortening of T| in Ramsey experiments. 

Room-temperature voltage noise was filtered using an anti-inductively wound coil of thin 
copper wire with a core of Eccosorb CRS-117 (~ 1 GHz cut-off), followed by two types of 
passive low-pass filters: 200 Hz second-order RC filters for DC biased lines, and 80 MHz 
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seventh-order Mini-Circuits LC filters for pulsed voltage lines (see supplementary section S4 
for measurements on the frequency response of the SET to a square wave applied to different 
gates.). The filter assemblies were placed in copper enclosures, filled with copper powder, 
and thermally anchored to the mixing chamber. DC voltages were applied using optoisolated 
and battery-powered voltage sources, connected to the cold filter box via twisted-pair wires. 
Voltage pulses were applied using an arbitrary waveform generator (LeCroy ArbStudio 1104), 
connected to the filter box via semi-rigid coaxial lines. ESR excitations were generated using 
an Agilent E8267D analog signal generator, and NMR excitations were produced by an Agilent 
MXG N5182A vector signal generator. Both excitation signals were combined using a power- 
combiner and fed to the MW antenna via a CuNi semi-rigid coaxial cable, with attenuators 
at the 1.5 K stage (10 dB) and the 20 mK stage (3 dB). The SET current was measured by a 
Eemto DLPCA-200 transimpedance amplifier at room temperature, followed by a floating-input 
voltage post-amplifier, a sixth-order low-pass Bessel filter, and acquired using a PCI digitiser 
card (AlazarTech ATS9440). 

S3 Data Acquisition Statistics 

Eor e“ experiments the state is always initialized spin-down and all of our plots were produced 
by taking the spin-up proportion from 100 — 200 single-shot measurement repetitions per point. 
For ^^P experiments, plots were produced by taking the nuclear flipping probability (no initial¬ 
ization to a certain state) from 41 measurement repetitions per point, and 50 electron readouts 
per nuclear spin readout. See Ref. (2) for more details on nuclear spin readout and control 
sequences. 
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S4 Frequency Response of the Electrostatic Gates 

All electrical gates are low-pass filtered to minimize the electron temperature. The top gate 
(TG), left barrier gate (LB) and the right barrier gate (RB) are filtered with 200 Hz second- 
order RC filters. The left donor “slow gate” (LDS) and the right donor “slow gate” (RDS) are 
filtered with nominally 10 kHz second-order RC filters. The left donor “fast gate” (LDF), the 
right donor “fast gate” (RDF), and the SET tuning gate (TGAC) are filtered with nominally 
80 MHz seventh-order Mini-Circuits LC filters (see Fig. S2 for schematic of the gate layout). 
However, the measured response of the SET to these different gates seems to differ from the en¬ 
gineered cutoffs. The reason for this is unknown, but could be explained by failed components 
or spurious conducting paths through the copper powder surrounding the lines inside the filter 
boxes. In Eig. SI we plot the current through the SET (colored regions indicate higher current) 
when a square wave is applied to the DS (LDS-i-RDS), DF (LDF-i-RDF) and the TGAC gates, 
respectively. For a frequency below the cutoff of the cable and filter, individual SET Coulomb 
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Figure SI: Frequency reponse. Frequency response of the SET current to a square wave 
applied to donor slow (DS), donor fast (DF) and the SET tuning gate (TGAC). The red solid 
lines are best estimates of the cutoff frequencies. 
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peaks are split into doublets. Above the cutoff frequencies the square wave is strongly attenu¬ 
ated and the Coulomb peaks merge back into single peaks. We estimate the cutoff frequencies 
of the different lines by comparing the measured data with theoretically modelled second-order 
RC filters, where Ri = 20 kfl, Ci = 1 pF, R 2 = 20 kfl, and C 2 = 1/(.R2/cutoff)- We find 
/cuiff = 0.6 MHz, = 3.6 MHz, and /^gac = j q mHz, and plot the corresponding 

frequency responses as red lines in Fig. SI. 

For the electrically-controlled measurements in the main text, a combination of DS, DF, and 
TGAC was used as “A-gate”. For the measurements in Figs. 2 & S7 the voltages applied to the 
gates were Vds = Va, ^df = Va, and IAgac = O.SV^. For the measurements in Figs. 3 & 4, 
and all other measurements in the supplementary information we reduced the voltage applied 
to the DS gate to improve the frequency response of the system. The applied voltages for these 
measurements were Vds = O.SV^, Vdf == Va, and V/gac = O.8I/4. 
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S5 Triangulation of Donor Position 


We can triangulate the position of the ion-implanted donor based on the techniques and methods 
introduced in Ref. (5). The triangulation is obtained by combining two different techniques, 
each one predicting a locus or several loci of donor locations compatible with a measurable 
physical property of the system. We use a classical, finite-element electrostatic simulation soft¬ 
ware (TCAD) (4), to model the electrostatic potentials in order to match the spin readout crite¬ 
rion (ground state energy of the donor-bound electron aligned with the Fermi level of the SET 
island), and a geometric capacitance extraction method (FASTCAP) (5), to match the measured 
capacitive coupling between the donor and the surrounding gate electrodes. 

For this device, we experimentally observed that the donor was more strongly capacitively 
coupled to TGAC than to the donor tuning gates LDF, LDF, RDS, and RDF. This indicates that 
the donor is positioned on the TGAC-side of the SET island. This may be due to misalignment 
of the implant window during fabrication, or to a donor belonging to the background doping 
of the epilayer. Matching the relative donor capacitances, obtained from charge stability ex¬ 
periments (see Table 1) with EASTCAP simulations, we obtain 4 loci for the possible donor 
location (see Figure S2). 

A 5th locus can be obtained from matching the spin readout criterion. The device had 
a threshold voltage of ~ 1.3 V. To faithfully describe the electrostatics of the experimental 
device, we include a negative interface charge density of Qox = —1-8 x 10^^ cm“^ in the TCAD 
model, necessary to match the threshold voltage (5). This charge density is consistent with 
estimates from deep level transient spectroscopic measurements (6). In a gated nanostructure, 
the donor can be susceptible to strain which can modify the conduction band energy {7, 8). 
Therefore, for our metrology, we choose a fairly large error bar for the spin readout criterion, 
accepting locations with = 45.6 ± 20 meV. The locus that matches the spin readout criterion 
is plotted in Fig. S2. We will see in Section S6 that the electric field at our final donor location 


S5 



Figure S2: Donor triangulation. Donor position extracted from a combination of donor ground 
state energy and relative donor capacitances to various gates. The combination of donor capac¬ 
itances and spin readout criterion aids to reduce the uncertainty of the donor position to ±4, 
±2.5 and ±3.5 nm in the three cartesian axis x, y and z. 


is ~ 6 MV/m, resulting in a slope of the conduction band of ~ 6 meV/nm. Hence, our large 
error bar in the spin readout criterion would only translate to a small error in donor position of 
~ 3 nm. 

The 5 shells in Fig. S2 intersect to within a region [50 ± 4, -31.5 ± 2.5, -8.5 ± 3.5], and 
represent the possible set of donor locations in the device. The donor is located under the right 
barrier gate RB, towards the TGAC gate. 
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Gates 

Relative Capacitance 

Cdonor—hB / Cdonor—TG AC 

0.65± 0.25 

{C'donor—LDS C!donor—KDS^/C'donor—TGAG 

1.53± 0.25 

{Cdonor—hDF Cdonor—KDF^/Cdonor—TGAG 

0.36± 0.075 

CTS Cdonor—island/Cdonor 

0.36± 0.06 


Table 1: Relative gate capacitances used for triangulation of the donor position. 

S6 Electric Field Simulations 

We can use the triangulated donor positions and our TCAD model of the Si/Si02 structure 
with the Al-gates (see section S5) to estimate the local electric field at the donor site. Fig. S3 
shows the magnitude and direction of the calculated electric field for = 50 mV (additional to 
Vlb = 0.92 V, Vrb = 0.92 V, Vtg = 1-798 V, Vos = 0.579 V, Vdf = 0.575 V, Vtgac = 0.45 V, 
compare also with section S4) inside a coordinate range that comprises the triangulated donor 
positions. Since the donor is located beneath the barrier gate, the positive bias applied to this 
gate (Vrb = 0.92 V) causes the electric field to point in negative z-direction, and the electron 


Electric Field for V = 50 mV (MV/m) 
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Figure S3: Electric field simulations. Calculated electric field for Va = 50 mV. 
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Electric Field Difference for A= 250 mV (MV/m) 
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Figure S4: Electric field simulations. Calculated electric field difference for AE = 


Ev^=S00mV — Evj^=b0mV- 


wavefunction is pulled towards the interface and away from the nucleus. This explains the 
very low initial hyperfine coupling A = 96.9 MHz, which is significantly different from the 
bulk value of 117.53 MHz. Increasing Vrb should pull the electron even further away from the 
nucleus reducing A even more. This means that we are starting in a situation slightly different 
from what is schematically depicted in Fig. IB, as the electron is already displaced towards the 
interface from the beginning. 

It may seem contradictory that, in experiment, we have measured a positive tuning parameter 
q;a = 0.91±0.07 MHz/V, meaning that a more positive Va leads to an increase in the hyperfine 
coupling. This is explained by the fact that the A-gate consists of the DS-gates, DF-gates, and 
the TGAC-gate (see section S4), which are all located some distance away from the donor. In 
this case an increase in Va will effectively compensate the electric field under the right barrier. 
We confirm this in Fig. S4, where we plot the change in electric field when Va is increased 


from Va = 0.50 mV to = 300 mV, i.e. AVa = 250 mV. The change in electric field is 
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Figure S5: Electric field simulations. Calculated electric field as a function of Va for possible 
donor locations. The change in electric field is very similar for all donor locations. 

clearly positive, reducing the magnitude of the vertical electric field when Va is increased, and 
therefore pushing the electron back towards the donor, resulting in a positive value of 

In Fig. S5 we plot the electric field E at possible donor locations as a function of Va- The 
stars correspond to the compensated plunge position (“Pulse ESR/NMR” position in Fig. 2A), 
that keeps the potential of the SET island constant with respect to the Fermi level of source and 
drain (9,1). The diamonds correspond to the electric fields calculated for the different values 
of Va- The absolute value of E varies significantly between the different donor positions, and 
converting Va into E would be subject to a large error. However, the tunability of E is very 
similar for all locations. Therefore, we fit this set of simulations to extract an average value for 
the tunability dE/dVA = —2.62(5) MVm“^A^, where the error is the standard deviation of the 
slopes of individual fits. This value is used to convert from Va to dE with good accuracy, and 
to calculate the “Electric Field Change”-axis of Fig. 2B. 
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S7 Atomistic Simulations of the Hyperfine Coupling 

The spin Hamiltonian of a donor electron spin S and nuclear spin I in an electrostatic 
potential (p and magnetic field Bq is given by: 

Hp-spin = 7e(0)S • Bo - 7nl ' Bq + ^(0)1 ' S (2) 

The first and second terms in Equation 2 are the electronic and nuclear Zeeman terms, the 
third term is the contact hyperfine interaction between the two spins, and 7e and 7 „ are the 
electron and nuclear gyromagnetic ratios, respectively. The contact hyperfine coupling (10) is 
expressed as >1(0) = ^ 7 e( 0 ) 7 „|' 0 (ro, (j))f, where |'0(ro, (j))f is the probability density of the 
electron wave function evaluated at the donor site tq. 

The relative tunability of the gyromagnetic ratio ay^/'ye(0) is two orders of magnitude 
smaller than the relative tunability of the hyperfine coupling aA/A(0) (refer to Fig. 2 and 
Ref. (11)). The hyperfine coupling relative to the bulk value (A(0) = 117.6 MHz) can, there¬ 
fore, be approximated as (3,12,13) 

A{<P) ^ |V^(ro,0)|^ 

>1(0) |0(ro,O)|^ 

We use a numeric implementation of tight binding - packaged as a software tool called 
NanoElectronic MOdeling-3D (NEMO-3D) (14,15) - to calculate the hyperfine coupling of 
the donor at the location triangulated in Section S5 and for the electric fields simulated in 
Section S6. Typical simulation domains of 30 nm x 30 nm x 30 nm consisting of approximately 
1.4 million atoms were considered. Each NEMO simulation with the above domain takes ~ 2 
hours, when run on a computing cluster with 48 processors. 

In Fig. S6A we plot the calculated hyperfine coupling for 1>4 = 50 mV for three different 
donor depths d at the triangulated position. Confinement by the interface enhances |0(ro, e)|^ 
beyond the bulk value for near shallow donors (d = 4.3 nm), resulting in a hyperfine coupling 
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greater than 117.6 MHz. However, for donors further away from the interface, the electric field 
pulls the electron away from the nucleus, resulting in a decrease in hyperfine coupling {13,16). 
The tunability of the hyperfine coupling A.A, i.e. how strongly A is modified when the electric 
field is changed, scales inversely proportional with the confinement of the electron. We plot 
A.A = ^v'^= 3 oomv “ ^yA= 50 mv fof three different donor depths in Fig. S6B. We notice that 
the calculated tunability of A is much larger than the experimentally measured one (compare 
Fig. 2). We believe this to be due to strain of the Si lattice, which also influences the hyperfine 
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Figure S6: Atomistic tight binding simulations of the hyperfine coupling for a donor at 
the location determined in Section S5, for the electric fields simulated in Section S6, and 
subject to lattice strain. A, B Hyperfine coupling and tunability of the hyperfine coupling as a 
function of donor depth under the Si02 in the unstrained lattice. C, D Hyperfine coupling and 
tunability of the hyperfine coupling as a function of lattice strain for a donor depth d = 4.3 nm 
under the Si02. The green lines indicate the experimental values and the green circles highlight 
simulations in good agreement with the experimental results. 
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coupling by distorting the donor wave function (i 7). 

We cannot quantify the strain in our device, but Thorbeck et al. (7) estimate that strain 
due to metal surface gates in nanoelectronic devices can be as much as ~ 0.1 %. Especially 
with the donor located directly under a barrier gate (see supplementary section S5) we expect 
significant strain at the donor location. We repeat the hyperfine calculations for the donor 
d = 4.3 nm below the interface for different values of homogeneous strain. Here, compressive 
strain (s < 0) in the xy-plane will lead to tensile strain (s > 0) in the z-direction. We plot 
the calculated hyperfine coupling for—0 .1% < s < 0.1% for strain applied to the xy-plane 
(blue squares) and the xz-plane (red circles) in Fig. S6C (Va = 50 mV). Any type of strain 
decreases the wavefunction overlap of electron and nucleus and reduces the hyperfine coupling 
A. In Fig. S6D, we plot the corresponding tunability of the hyperfine coupling AA (for AVa = 
250 mV). For certain values of s both the absolute value of A and the tunability A A are in 
good agreement with the experimental values of A^^p = 96.9 MHz and AA^xp = 0.228 MHz 
(highlighted with green circles). While the presence of strain in the device seems to be able to 
match the simulated values with the experimental values, we cannot conclusively attribute the 
observed effects to a specific value of strain. Further and more detailed simulations would be 
required to untangle the effects of electric field, strain and the Si 02 interface on the hyperfine 
coupling. 

Electric fields, strain and the Si02 interface have profound consequences for future multi¬ 
qubit devices as they can lead to a strong variability in the hyperfine coupling A and, therefore, 
to distinctively different qubit resonance frequencies. It can, however, be expected that multi¬ 
qubit devices with deterministically positioned donors will have a smaller variability in the 
hyperfine coupling as the donors can be placed at locations with little or at least similar strain. 
Furthermore, a gate layout optimized for maximizing the Stark shift should give a tunability 
of the resonance frequencies of a few MHz and be sufficient for the operation of a multi-qubit 
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quantum computer with a monochromatic global microwave field. In a scenario where the 
variability is too large or the tuning range too small, the global microwave could be operated as 
multi-tone continuous-wave driving field at regular frequency spacing. 
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S8 Electron Spin Resonance Spectrum & Rabi Oscillations - Experiment 
& Theory 

In Fig. S7B we show a gate-controlled measurement of the ESR spectrum, obtained by shifting 
Pe with Va- The nuclear spin is in the Ifl) state for the duration of the experiment, and we apply 
a continuous-wave MW driving field at frequency i^mw = “ 200 kHz. Here, PI 2 is the ESR 

frequency for jfi^) obtained from a conventional pulsed Ramsey experiment at the compensated 
plunge position (“Pulse ESR/NMR” position in Fig. 2A). We start the sequence (see Fig. S7A) 
at 1/4 = 0 V to load an electron in the |4.) state by spin-dependent tunneling (9) (“Read/Init.” 
position in Fig. 2A). We then apply a positive Va pulse to shift Pe 2 {VA) towards pmw for the 
duration tp, which results in a coherent manipulation of the spin. We then pulse Va back to 0 V 
to perform single-shot readout of the electron spin (9), and we repeat the whole sequence 200 
times to extract the electron spin-up fraction P-^{Va). For Va = Vr ~ 155 mV, an increased 
count of spin |t) electrons indicates that Pe 2 {VA) becomes resonant with the MW source. By 
optimizing the duration of the gate pulse we can ensure that the electron spin undergoes a tt- 
rotation while on resonance, yielding ~ 1. 

The maximum P^ occurs at a frequency shift APe = —217 kHz instead of the expected 
—200 kHz. This shift and the presence of side lobes at larger Va, is caused by the finite time re¬ 
sponse of the electrical control gates, which are low-pass filtered to minimize electron heating 
(see supplementary section S4 for details). When measuring the gate-controlled ESR spec¬ 
trum beyond pmw (Va > Vr), we sweep Pei(VA) through the resonance and back again. The 
limited bandwidth of the gates causes slow crossings through the resonance condition, result¬ 
ing effectively in a Landau-Zener-Stiickelberg interferometry experiment (18). Time-evolution 
simulations of the whole sequence (red line in Fig. S7B), taking into account the bandwidths 
of the different gates, show excellent agreement with the measured data. In a setup optimized 
for electrical control, the rise-time of the A-gate should be chosen much shorter than the qubit 
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Figure S7: Electrically-controlled ESR spectrum. A, Schematic of the sequence used to 
measure the electrically-controlled ESR spectrum. B, ESR spectrum obtained by electrically 
tuning the ESR transition (for Ifl)) in resonance with the continuous-wave MW source (vmw = 
Pe 2 — 200 kHz) with a short voltage pulse. Time evolution simulations (red line) are in excellent 
agreement with the data. 


Rabi period while on resonance. 

In the following paragraphs we present the theoretical model that we have set up to describe 
the temporal response of the donor system to a voltage pulse on the A-gate, and the temporal 
evolution of spin qubits when they are tuned into resonance with a CW magnetic driving field. 
We model the experimental data using the density matrix formalism with the Hamiltonian 

H{t) = ^hAp{t)a^ + ^hfloa^, (4) 

where Av{t) = Ve2 { t ) — i^mw is the detuning between the ESR transition and the MW source, 
and Qo = 23.8 kHz is the Rabi frequency for a Bi = 0.85 /uT, which gives a 7r-pulse length of 
21 /iS. The output power of the MW source Pmw = “22 dBm was chosen very low for these 
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experiments to reduce the power broadening of the ESR line. The dephasing time T 2 = 970 /US 
of the electron spin (see Fig. 3H) is included in the master equation of the Lindblad form (19) 

| = (5) 

where 


^(p) = 7^ (2cr^p(T^ - cr^cr^p - pa^a^) 

Z1 2 

2 

= -T^<7xP- (6) 

J-2 

We then use the equation of motion (5) to numerically compute the time evolution of an electron 
initialized in |4.). 

From the measurements in Fig. 2B, we know dPe 2 /dVA and can allocate the total shift in Ve 2 
to DS, DF, and TGAC in accordance to their respective relative capacitive couplings to the donor 
(see section S5) and temporally filtered to their respective cutoff frequencies (see section S4). 
This allows us to generate the detuning trace Az/(t) for any measurement and detuning sequence 
that we want to model: 


Au{t) 


{dv,2/dVA) 


Ds/C^donor—TGAC “1“ Dr/C^cionor— TGAC + 1 

{iCdoncn — Ds/C'<ionor-TGAc)G^^(f) + (C'donor-DF/C'donor-TGAc)G^^(f) + Vj 

lA+McTl + 0.36Kf"(i) + . 


TGAC 


(7) 

(*)) 


Here, and are the low-pass filtered (/cutoff = O-® MHz, = 

3.6 MHz, and — 7.0 MHz) voltage traces that are applied to the gates. 

Fig. S8A shows Av{t) for an electrically-controlled ESR measurement. The time evolution 
simulation starts with an electron in the |/) state at Av = 274 kHz for the first 10 ps. Av 
is then changed to —25 kHz (in this specific example) for a time of tpuise = 21 ps (yr-pulse) 
before it is tuned back to 274 kHz for 16 ps. Fig. S8B,C,D show the time evolution of the 
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Figure S8: Time evolution simulations of the electrically controlled ESR spectrum. Simu¬ 
lations depicting the time evolution of the electron spin around the Bloch sphere for a detuning 
pulse of tpuise = 21 fxs and amplitude 249 kHz that tunes the ESR transition to 25 kHz lower 
frequency than the MW drive. 


electron spin during this sequence. The final z-orientation of the electron spin is then plotted 
in Fig. S8E, and the whole calculation is repeated for different detuning values to build up 
the entire ESR spectrum (comp. Fig S8E and Fig. S7B). The calculated spectrum is in very 
good agreement with the experimental data when taking into account non-unity readout fidelity 
(0.85 %) and non-zero background counts (0.12). The sidelobes at Ai/ < 0 are well reproduced. 
They appear when Ue 2 is tuned through umw at the beginning of the pulse and back at the end. 
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The limited bandwidth of the A-gate causes slow crossings through the resonance condition, 
resulting effectively in a Landau-Zener-Stiickelberg (LZS) interferometry experiment (18). 

The measurements of Fig. S7B and the calculations of Fig. S8 can be also be performed as an 
electrically-controlled Rabi experiment, where the final electron spin orientation is measured as 
a function of the length of the voltage pulse. In Fig. S9 we plot the experimentally determined, 
electrically-controlled Rabi spectrum in the left panel, and the simulated one in the right panel. 
Again, the simulations are in excellent agreement with the experimental data and the evolution 
of the LZS sidelobes is well reproduced. 


Experiment Time Evolution Simulations 



Shift in ESR Frequency (kHz) Shift in ESR Frequency (kHz) 


Figure S9: Electrically-controlled Rabi spectrum. Experimental data and time evolution 
simulations on the electrically controlled Rabi spectrum. 
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S9 Ramsey Experiments 

The ability to apply VA{t) sequences to the A-gate allows us to perform arbitrary qubit control 
sequences as already demonstrated in Fig. 3 and Fig. 4. This supplies us with an alternative way 
to measure the frequency shift induced by Va- While the measurement in Fig. 2B was conducted 
changing Va during the control phase and pulsing the MW source, we can also conduct an 
electrically-controlled Ramsey experiment as introduced in Fig. 3D. In Fig. SlOA we plot a 
series of Ramsey measurements, where AVa = i.e. the detuning of the spin 

transition during the wait time, was changed from one measurement to the next. The frequency 
of the Ramsey oscillations corresponds to the detuning, and we plot the extracted values as 
function of AV^ in Fig. SIOB. The slope of the linear fit is dPe 2 /dAVA = —0.80±0.04 MHzW. 

Fig. Sll shows a similar data set for the neutral nucleus. The slope of the linear fit is 
dPni/dAVA = 0.30 ± 0.01 MHz/V, which allows us to calculate the tuning parameters a a = 
dA/dAVA = 2dPni/dAVA = 0.60±0.02 MHzW, and a^^Bo = d'^^Bo/dAVA = dUe 2 /dAVA — 
dPni/dAVA = -1.10 ± 0.05 MHzW, with = d-f^jdAVA = -0.71 ± 0.04 MHzW/T at 
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Figure SIO: Electrically controlled electron Ramseys. A Set of EC Ramsey experiments 
for different voltage pulses AVa conducted on the e“. B Corresponding frequency shift as a 
function of AlA- 
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Figure S11: Electrically controlled nuclear Ramseys. A Set of EC Ramsey experiments for 
different voltage pulses AVa conducted on the neutral nucleus. B Corresponding frequency 
shift as a function of AVa- 


Bq ^ 1.55 T. All the values extracted in this section are slightly smaller than the values extracted 
from Fig. 2. This is because we reduced the voltage applied to the DS gate (V^s = 0.5Va 
compared to V^s = Va for the measurements in Fig. 2) to improve the frequency response of 
the system for all electrically-controlled measurements involving pulsing sequences (see also 
section S4). 
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SIO Coherence Times Measurements 


In addition to the Hahn echo and CPMG measurements introduced and presented in Fig. 3G,J,K, 
we performed Ramsey experiments to extract T 2 and dynamical decoupling experiments with 
a different number of refocussing pulses. The result of these measurements is presented in 
Fig. S12 for the e“ and in Fig. S13 for the neutral nucleus. Furthermore, the measurements 
on the e“ have been performed as both electrically controlled and conventional, pulsed spin 
resonance experiments to directly compare these two measurements methods. The Ramsey 
experiments are displayed in Fig. S12A,B, the measured CPMG decay traces and their fits are 


Electrically Controlled Spin Resonance Experiments 
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Conventional Pulsed Spin Resonance Experiments 



Figure SI2: Electrically controlled electron coherence times. A,B Ramsey experiments to 
extract T 2 for the electron qubit. C,D CPMG dynamical decoupling decay traces for different 
numbers of refocussing pulses. E,r Extracted coherence times as a function of number 

of refocussing pulses N. 
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Figure S13: Electrically controlled nuclear coherence times. A Set of EC CPMG dynamical 
decoupling performed on the neutral nucleus. B Extracted coherence times as a 

function of number of decoupling pulses. 


displayed in Fig. S12C,D and Fig. S13A, and the extracted decay times are plotted as a function 
of number of CPMG pulses in Fig. S12E,F and Fig. S13B. 

The electron spin coherence times obtained for electrically controlled and pulsed experi¬ 
ments are identical within the error bars, indicating that there is no additional source of de¬ 
coherence introduced by performing electrically-controlled measurements. The free induction 
decay time is T 2 = 220fxs (see Fig. S12A,B) and the extended coherence times reach values 
of = 10 ms for 32 decoupling pulses (see Fig. S12C-F) for both experimental meth¬ 

ods. The electron spin coherence times can be fitted with yCPMC js^a/ia+i) _ 7 ^ 0 . 67 ^ where 
a = 2.0 gives information about a colored noise spectrum with S{lj) oc l/a;^ ° in good agree¬ 
ment with the data obtained by noise spectroscopy with pulsed-MW control sequences with 
j,puiseCPMG ^ o.93A^° (a = 2.3). This is another strong indicator, that the electrically- 
controlled measurements do not introduce significant additional noise. 

The coherence times of the neutral nucleus can be fitted with oc = jyo 56 ^ 

with a = 1.3 that hints towards l//-type noise. In contrast to the measurements presented 
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here, in measurements with pulse-RF control sequences we were not able to significantly extend 
with dynamical decoupling sequences (see Ref. (20)). We believe that this is caused by 
a sensitivity of the donor system to the RF radiation, possibly by a process as simple as heating. 
For the pulse-RF technique, the RF power was chosen higher and pulsing of the RF source will 
lead to a non-equilibrium state. On the other hand for the electrically-controlled technique, the 
RF power was chosen lower (to reduce power broadening of the NMR transition) and the CW 
RF drive was leading to a steady-state situation. 
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